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Abstract. The density in the core of neutron stars can reach values of about 5 to 10 times nuclear matter 
saturation density. It is, therefore, a natural assumption that hadrons may have dissolved into quarks under 
such conditions, forming a hybrid star. This star will have an outer region of hadronic matter and a core of 
quark matter or even a mixed state of hadrons and quarks. In order to investigate such phases, we discuss 
different model approaches that can be used in the study of compact stars as well as being applicable to 
a wider range of temperatures and densities. One major model ingredient, the role of quark interactions 
in the stability of massive hybrid stars is discussed. In this context, possible conflicts with lattice QCD 
simulations are investigated. 

PACS. 26.60.-c Nuclear matter aspects of neutron stars - 26.60.Kp Equations of state of neutron-star 
matter - 25.75.Nq Quark deconfinement, quark-gluon plasma production, and phase transitions 


1 Introduction 


The physics of neutron stars is one of the central areas 
of research in nuclear astrophysics. Neutron stars rep¬ 
resent a unique environment for investigating extremely 
dense strongly interacting matter at relatively low tem¬ 
peratures. This renders neutron star physics as an ideal 
complement to the efforts using relativistic heavy-ion col¬ 
lisions in the study of matter under extreme conditions. In 
the latter case the fireball created in the collision zone has 
a very high temperature, thus sampling a different regime 
of QCD matter than neutron stars. 

Using equation of states of purely hadronic matter, 
typical maximum densities in the center of the stars range 
from about 5 to 10 times nuclear matter saturation den¬ 
sity. As these numbers correspond to close packing from 
a purely geometrical point of view, it is very natural to 
assume that at this point (and perhaps much earlier) the 
baryons have dissolved into their quark components. Such 
an object is termed hybrid star, i. e. a star with a hadronic 
outer core surrounding a quark, or mixed hadron-quark, 
inner region ^[2||3IBl5ll6ll3[^E| riO[ITT] . Depending on the 
equation of state used to describe the quark phase, very 
different masses and radii are obtained for stars when com¬ 
pared to a purely baryonic star. We will discuss this issue 
in detail in the following sections. 


2 Modeling Compact Stars 

In order to investigate the properties of hybrid stars, we 
first introduce the formulation of the hadronic chiral mean 
field (CMF) model. It is based on a SU(3) extension of a 
sigma-omega chiral model in a non-linear realization of 
chiral symmetry, including all lowest hadronic multiplets 
[HE]. The advantage of such an approach is that it also 
incorporates dynamic mass generation, which allows for a 
realistic description of chiral symmetry restoration, while 
being in very good agreement of nuclear and astrophysics 
properties. 


2.1 The Hadronic Model 

The CMF model includes nucleons, hyperons as well as 
non-strange and strange mesons. As in the usual relativis¬ 
tic mean field approach, the baryons interact via meson 
fields. In addition, in the CMF approach also the bary¬ 
onic masses are largely generated by the interaction with 
the scalar fields, except for a small explicit mass term. 
The specific form of the potential of the scalar fields leads 
to non-vanishing vacuum expectation values, thereby pro¬ 
ducing a dynamic baryon mass via the baryon-meson cou¬ 
pling, analogously to chiral sigma models. 

The terms in the Lagrangian density relevant for the 
following discussions read [HE! 
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^Kin contains the kinetic energy of the fields. The interac¬ 
tion term £i nt defines the linear interaction between me¬ 
son fields and baryons. It reads explicitly 

Ant = - X] rfihoigiuV + 9i4>(t> + QipTzp) + M*]lpi . (2) 

i 

The baryons i interact with the vector mesons u (isoscalar), 
p (isovector), and (j) (isoscalar with hidden strangeness). 
The various coupling constants qbm are fitted as discussed 
below based on SU(6) symmetry [12]. The effective masses 
M* are generated by the coupling of the baryons to the 
scalar mesons, reading 

M* = gi a (7 + gisrsS + + Mq. . (3) 

These expressions include couplings to the various scalar 
fields, the isosclar cr, the isovector S and the field with hid¬ 
den strangeness £. The couplings are connected via SU(3) 
relations m- 

The scalar fields have self-interaction terms that are 
responsible for generating non-vanishing vacuum expecta¬ 
tion values for the a and ( that, in quark language, cor¬ 
respond to the scalar vacuum condensates of non-strange 
and strange quarks, respectively. Using SU(3) invariant 
terms for the interactions, the Lagrangian term is given 

by 


£seif = ~^(m 2 u> 2 + m 2 p p 2 + m \<\> 2 ) 

+ k 0 (a 2 + C 2 + S 2 ) + ki (cr 2 + C 2 + S 2 ) 2 
+ + 3a 2 S 2 + C 4 ) + k^(a 2 — <5 2 )C 


+ Aq 

+ #4 


In 


(cr 2 - (5 2 )C 

°o(o 


to 4 + 
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+ 3u 2 (f> 2 + 


4w 3 4> 


+ 



(4) 


Finally, an explicit symmetry breaking term that also gen¬ 
erates masses for the pseudoscalar mesons is included 


£sb = + (\/2 m 2 k f k - C • (5) 

The equations of motion for hadronic matter at zero or 
finite temperature are determined from extremizing the 
grand canonical potential 


&/V = - Li n t ~ Lself ~ Lsb ~ Ly ac 
T^E, J2W fo Fi d3k ln (l ± e-*W(‘W)), (6) 

including a heat bath of quasiparticles. This model ap¬ 
proach has been used in many applications in relativistic 
heavy-ion physics mm, nuclear structure calculations 
fT6lfl7l . as well as astrophysical studies [T81H9] . Fig. [l] 
shows a sample of these equations of state at vanishing 
temperature m- In-between the equation of state of nu¬ 
cleonic matter and the relatively soft equation of state for 



Fig. 1. Equation of state of the CMF hadronic model at T = 
0. Each curve shows the equation of state for different cases: 
isospin symmetric nuclear matter and star matter assuming 
different degrees of freedom. 


isospin symmetric matter, the figure also incorporates re¬ 
sults taking into account the hyperon states as well as an 
extension of the calculation that also includes the spin-3/2 
decuplet of baryon resonances. As can be observed, adding 
further degrees of freedom leads to a softening of the equa¬ 
tion of state [T8lfT9] . Note that for the isospin symmetric 
nuclear matter calculation net strangeness is set to zero, 
whereas in general stellar matter is calculated assuming 
beta equilibrium and charge neutrality. 


2.2 Hybrid stars 

As has been argued in the introduction, it is very natural 
to assume that, at some density, quark degrees of freedom 
appear inside compact stars. Before discussing the exten¬ 
sion of the hadronic model to include quarks, let us look 
at some general aspects of hybrid star solutions. For this 
purpose, we combine an often-used hadronic model (G300 
[20] ) with a simple MIT bag quark equation of state EU 
[22] . In this approach we can adjust the quark phase bag 
pressure B. The choice of B essentially determines at what 
chemical potential the phase transition to quark matter 
occurs and, connected to this, it determines the strength 
of the first-order phase transition. 

Fig. m shows that for larger values of the bag pres¬ 
sure the transition point is shifted to higher densities and 
the jump in the energy density also significantly increases. 
Both of these have direct consequences for hybrid star 
solutions. The curves in Fig. [3] are the solutions of the 
Tolman-Oppenheimer-Volkov (TOY) equations for static 
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Fig. 2. Equation of state of the G300/MIT hybrid model. 
The various curves correspond to different choices of the vac¬ 
uum pressure with values R 1//4 = 160,165,167.5,170,180 MeV 
(larger values generating later transitions). 



Fig. 3. TOV solutions for the equations of state shown in Fig. 
[2] Strong first-order phase transitions with large latent energies 
lead to stellar instability when quarks appear. 


spherically symmetric stars mm for this set of equa¬ 
tions of state. The kinks in the curves in Fig. [2] mark the 
onsets of the quark phase. The general effect of this insta¬ 
bility is a reduction of the maximum star masses as well 
as a shift of the radius of these stars to smaller values. 
In the case of very strong transitions, the onset can even 
lead to a removal of all stars with a quark core from the 
range of stable solutions. This is the case for the equation 
of state generated with a bag pressure B = (180MeV) 4 . 
For some cases, also two maxima develop, which would 
correspond to twin star solutions, meaning two branches 
of stable hadronic and hybrid stars, respectively I25ll26l 

mmm- 

This general behaviour can be modified by including 
interactions between the quarks as shown in Figj4j Fol¬ 
lowing the simplified procedure outlined in Ref. [30], the 



Fig. 4. As Fig. [3] but effectively including a repulsive inter¬ 
action to the quark equation of state. Results for coupling 
strengths a = 0.1,0.3,0.7 are shown. Large repulsion leads 
to larger star masses and weaker phase transitions. 

modified quark phase is described by the grand canonical 
potential 

Gq = + B • ( 7 ) 

Here, the additional term proportional to the strong inter¬ 
action coupling strengh a s contains the repulsive interac¬ 
tion between the quarks. As was pointed out in Ref. m, 
adjusting the interaction strength accordingly, the tran¬ 
sition between quarks and hadrons becomes weak, effec¬ 
tively leading to similar masses and radii for the hadronic 
and hybrid stars. In this way, a strong repulsive quark in¬ 
teraction helps to stabilize hybrid stars. A general discus¬ 
sion of the stability of hybrid stars, using model-independent 
parameterized equations of state, can be found in Ref. m- 
However, this procedure of stabilizing hybrid stars creates 
different problems, as it will be briefly discussed in Section 

m 

2.3 Unifying hadron and quark models 

As the model approach under discussion is supposed to 
be (and has been) applied to a large range of temper¬ 
atures and densities, including conditions that exist in 
the fireball of heavy-ion collisions, its general phase struc¬ 
ture should be in accordance with the known features of 
the QCD phase diagram. The well-established properties 
are (i) bound nuclear matter, i.e. a first-order liquid-gas 
phase transition at the nuclear matter saturation density 
of about 0.15 fm -3 and (ii) a smooth crossover for vanish¬ 
ing net baryon density at a temperature of about 160 MeV 
as it has been unambiguously established in lattice QCD 
simulations [321133] . Thus, a simple connection of two sep¬ 
arate models of the hadronic and the quark phases, as it 
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Fig. 5. Phase diagram of the quark-extended CMF model. 
Apart from the nuclear liquid-gas first-order phase transition, 
the chiral/deconfinement transition exhibits a critical end¬ 
point at small chemical potentials. Results for isospin sym¬ 
metric and star matter are shown. 

was used in the previous section, is not sufficient to ful¬ 
fil these requirements. Connecting both equations of state 
will necessarily lead to a first-order transition for all condi¬ 
tions (except perhaps for a single point, i.e. chemical po¬ 
tential, by fine-tuning parameters). Therefore, a unified 
model including hadronic and quark degrees of freedom 
is required. In [5] such an approach has been developed 
for the first time (see also Refs. |34ll35] . In this formula¬ 
tion, quarks and baryons interact via the scalar and vector 
meson condensates. In addition, in analogy to the quark 
Polyakov loop - NJL (PNJL) model, an effective field <P 
is introduced that effectively describes the order param¬ 
eter for the deconfinement transition [36,37 . This field 
has a potential energy term, which is fitted in order to re¬ 
produce the thermodynamic and Polyakov loop measure¬ 
ments in lattice QCD simulations. A simplified potential 
of this type was proposed in Ref. [5] 

= (a 0 T 4 + ai/i 4 + a 2 T 2 /i 2 )^ 2 

+a 3 T 4 log (1 - 6<Z> 2 + 8<Z> 3 - 3<£ 4 ). (8) 

The constants ao and a 3 were fitted to obtain a reasonable 
agreement with lattice results at zero chemical potential. 
The potential term U also contains terms that depend on 
the chemical potential. They can be adjusted to gener¬ 
ate a first-order phase transition at zero temperature to 
a deconfined state at a chosen baryon chemical potential. 
Further, they can be adjusted to lead to a phase diagram 
with a critical end-point of the first-order phase transition 
line that is in agreement with a particular lattice QCD 
calculation [38]. The resulting phase diagram of such an 
approach is shown in Fig. [5] The figure exhibits a transi¬ 
tion line with the end-point at small chemical potential. 
It also shows the nuclear liquid-gas phase transition. The 
two slightly shifted lines mark the difference of the tran¬ 
sition for isospin-symmetric and star matter. 

Another ingredient of this model approach is the cou¬ 
pling of the <3 field to the masses of the baryons and 



Fig. 6. Equation of state for star matter in the quark-extended 
CMF model at zero temperature. Two assumptions, local 
(Maxwell) and global (Gibbs) charge neutrality, are consid¬ 
ered. 

quarks. This is a simple phenomenological way of ensuring 
that, for small value of the field, quarks are not populated 
and vice versa in the case of baryons. The effective masses 
read explicitly (for details, see Ref. [5]): 

Mg = gBaC + gBs^S + gBcC + M 0b + gB<p@ 2 , (9) 

M* = g q(J a + g q 8Ts5 + g q ^( + M 0q + g q <p( 1 — 3). 

( 10 ) 

With this model in hand, it is now possible to calcu¬ 
late hybrid stars as well as quark-hadron matter at high 
temperatures and/or densities in general. Using the model 
parameters as given in Ref. [51139] , the equation of state of 
star matter can be determined. The result for vanishing 
temperature is shown in Fig. [6| As can be seen in this 
figure, the equation of state has a strong first-order phase 
transition. The figure includes two different treatments of 
the transition: first, local charge neutrality in both phases 
is assumed, which corresponds to a Maxwell construction, 
and a transition that happens at a fixed pressure. If one 
relaxes this condition to a global charge neutrality, with 
opposite charges in quark and hadronic phases, one ob¬ 
tains a range with a mixed phase and varying pressure. 
Which of these scenarios is correct, depends on the QCD 
surface tension which has to be taken into account in a 
mixed phase with bubbles of different phases. Unfortu¬ 
nately the value of the surface tension, not to mention its 
variation with density, is largely uncertain mm- 

The equations of state from Fig. [6] can then be used 
to solve the TOV equations. The obtained results are pre¬ 
sented in Fig. [7] The full line in the large figure is the 
result for stars without taking into account quarks but 
including hyperons. The main reason that the hadronic 
star including hyperons (hyper star) still has a large max¬ 
imum mass, which is only reduced by about O.O6M 0 [T5] 
02], arises from the fact that in this full ffavor-SU(3) ap¬ 
proach the hyperons experience a strong repulsive force 
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Fig. 7. Mass-radius diagram of the quark-hadron model from 
Section |2.3| The strong first-order phase transition leads to 
unstable hybrid stars using Maxwell construction. With the 
Gibbs procedure, hybrid stars with a core of a mixed state of 
quarks and baryons become stable. 


from <p and uo meson exchange. As an illustration of this 
feature, Fig.[8]shows the resulting maximum mass and the 
strangeness content f s (defined as the average strangeness 
per baryon) in the center of the star as function of the 
scaled hyperon-omega meson interaction (note that one 
could alternatively also scale the interaction with the <p 
meson). Decreasing artificially the coupling shows a sub¬ 
stantial increase in strange particles and the accompany¬ 
ing drop in maximum mass. 

Going back to the discussion of figure [7| including the 
quarks in the calculation, the resulting hybrid star so¬ 
lution features a kink at around 2 solar masses. As we 
discussed in more general terms in Section |2.2[ the latent 
energy in the phase transition adds a large amount of en¬ 
ergy density to the star matter when quark states begin to 
be populated, but not much more pressure. Therefore, us¬ 
ing the Maxwell-constructed equation of state, the hybrid 
star becomes unstable with respect to collapse. On the 
other hand, using the Gibbs construction, a mixed phase 
of quarks and baryons survives in the core of the star as is 
indicated in the inset of the figure. This region can extend 
up to 2 km in the center of the star. 

It is interesting to note that, although the maximum 
central density shifts from about 6po to 4po (in the Maxwell 
case) due to the instability introduced by the quark phase 
there is only a moderate reduction in maximum mass of 
about 0.1 solar masses. The reason for this behaviour is 
illustrated in Fig. [9] Here, the star masses are plotted as 
function of central density for the purely hadronic stars. 
The figure shows that the mass is quite insensitive to the 




Fig. 8. Stellar mass and central strangeness fraction f s as func¬ 
tion of the hyperon-vector meson interaction strength. Reduc¬ 
ing the coupling substantially increases the strangeness frac¬ 
tion and reduces star masses. 

central density for a wide range of densities below the 
maximum mass value. Therefore, in this case even a re¬ 
duction of the value of the central density by about 30 
percent does not appreciably change the maximum star 
mass. 


2.4 Excluded volume corrections 

In the development of the quark-hadron approaches as dis¬ 
cussed in this article, a second ansatz that incorporates a 
different mechanism in transiting from baryons to quarks 
at high densities and temperatures was developed. The 
hadronic model used is the same as the one discussed in 
Section |2.1| In the calculation presented here, the expres¬ 
sion for the Polyakov loop potential is the one used in a 
number of PNJL calculations 36j. 37 : 

U = -fa(T)<2><r 

+ b(T)ln[ 1 - 6<2><T + 4(<2> 3 <r 3 ) - 3(<2><T ) 2 ] . (11) 

The particular logarithmic term, which is also present in 
Eq. ©> corresponds to the Haar measure of integrating 
over the SU(3) Polyakov loop matrix. The temperature- 
dependent couplings are defined as a(T ) = aoT 4 +aiT 0 T 3 + 
g^XqT 2 and b(T) = b^T^T. The numerical values of the 
various parameters are adopted from Ref. [36] . This choice 
ensures a correct first-order phase transition in the pure 
gauge sector, i.e. for the theory including only the <P field. 
Again, following the spirit of the PNJL model, the distri¬ 
bution functions of the quarks and antiquarks contain the 
































6 


S. Schramm et ah: Modelling Hybrid Stars in Quark-Hadron Approaches 



Fig. 9. Star masses as function of their central density p c for 
the hadronic CMF model. As can be seen for masses around 
2 solar masses, the actual star mass is relatively insensitive to 
p c in this region. 

Polyakov loop field 

Q o = - T Y, (2^3 / rf3fcln i 1 + ^ ex P - T ’ 

( 12 ) 

and 

jeQ v ' J v 7 

(13) 

where 7 \ stands for the spin and color degrees of freedom. 
This ensures suppression of quarks in the confined state 
(@ = 0). In the case of hadrons, an excluded volume term 
is introduced that effectively mimics short-range repulsion 
of overlapping hadrons. Here, we follow a very simple ap¬ 
proach, assuming different excluded volumes for quarks, 
baryons, and mesons of the form 

VQuark — 0 i 
VBaryon ^ 5 

VMeson — V/CL , 

(14) 

with a value v = 0.64 fm 3 . The parameter a takes into ac¬ 
count that the mesonic volume is smaller than the one for 
baryons. We assume it to be a = 8 , which implies a meson 
radius of half the size of the baryonic one. In this basic 
implementation, we neglect possible Lorentz contraction 
effects for the excluded volumes, which were discussed in 
Ref. [431144] . or density or temperature-dependent hadron 


sizes mm- For thermodynamic consistency, the modi¬ 
fied chemical potential Jli of particle species i is given by 

Mi = Hi - Vi P , (15) 

where P is the sum over all partial pressures. In this for¬ 
malism the thermodynamic quantities are calculated with 
respect to the temperature T and the modified chemical 
potentials Jli. Thermodynamic consistency is fulfilled by 
multiplying the densities (energy density number den¬ 
sity pi , entropy density si), calculated with the values 
by a factor /, which is the ratio of the total volume V to 
the unoccupied reduced volume V' 

/ = y = (i + XA^) -1 • ( 16 ) 

i 

Thus, for instance, the energy density then reads 

e = Ylf % ■ ( 17 ) 

i 

Some results of this model for vanishing baryon chemical 
potential are collected in Fig. [l0| The plots compare model 
results to lattice QCD calculations [47]. The left panel 
shows the temperature dependence of the Polyakov loop 
and the scalar condensate, whereas the right panel depicts 
the densities of quarks and hadrons (summing particle and 
anti-particle densities). As the transition in this regime is 
a smooth crossover, it is rather natural to assume that 
there is a mixed state, and not a mixture of two separate 
phases, of hadrons and quarks as it is observed in this 
model. 

When calculating the phase structure of isospin sym¬ 
metric matter for a wide range of temperatures and densi¬ 
ties, this model does not lead to a first-order phase transi¬ 
tion (in addition to the liquid-gas phase transition). This 
feature does not violate any established constraints, as the 
type of phase transition of matter at high density is un¬ 
known. However, in the case of charge-neutral matter in 
beta equilibrium, the situation is different. More precisely, 
some solutions for different interaction strengths of the 
non-strange quarks to the uj field, g quj , and of the strange 
quark to the corresponding (j) field, g s (with the relative 
strength £ = g S( f)/g quJ ) produce a first-order transition to a 
strangeness-enriched phase for star matter. This is shown 
in Fig.[lT]for the strange scalar field £. The corresponding 
solutions for the star families shown in Fig. [12] demon¬ 
strate that such a transition might lead to a small range 
of second family “twin” stars (see also Ref. [29] for the 
case of very heavy twin solutions). Looking at the differ¬ 
ent particle densities as is shown in Fig. [13] one can see 
a steep rise in the s quark population beyond the phase 
transition. Note that even stars that do not reach the crit¬ 
ical density of the phase transition already contain quarks, 
being also hybrid stars. 

Twin star solutions have recently been revisited in the 
literature as they pose a particular solution for the so- 
called “hyperon puzzle”. As the appearance of quark mat¬ 
ter suppresses the amount of hyperons at high density, this 
scenario might allow for small, but massive stars. 









S. Schramm et ah: Modelling Hybrid Stars in Quark-Hadron Approaches 


7 




temperature [MeV] temperature [MeV] 

Fig. 10. Model results including excluded volume corrections at vanishing chemical potential. The figure contains plots of the 
temperature dependence of the effective Polyakov loop field the scalar field cr, and the densities of quarks and hadrons with 
data points from lattice calculations El- 



Fig. 11. Scalar field £ with hidden strangeness normalized by 
its vacuum value. Depending on the quark interaction strength 
£ a first-order transition might occur for star matter. 



Fig. 12. TOV solutions for hybrid stars using different 
strengths of the quark-ca and quark-0 couplings. For specific 
choices, a twin-star solution might occur. 


One very interesting and promising venue in the study 
of compact stars is the investigation of their thermal evo¬ 
lution. This evolution is strongly dependent on the particle 
composition of the star, which makes it an ideal tool to 
gather information on the actual particle content of the 
star. The equations that describe the stellar thermal evo¬ 
lution in a general relativistic setting read [ 48] 


<9(Ze 2 ^) 

dm 


1 

pyjl — 2 m/r 



~hc, 


d(Te^) 




d(Te*) _ (le*) 

dm 167r 2 r 4 ftpy/l — 2 m/r 

Here, the variables r, p{r) and m(r ) denote the radial 
distance from the center of the star, the energy density, 
and the stellar mass, respectively. The various thermal 
properties that enter the evolution are the temperature 
distribution T(r, £), luminosity Z(r, t), the neutrino emis- 
sivity e v (r, T), thermal conductivity ft(r, T) as well as the 
specific heat c v (r,T). The equations have to be supple¬ 


st 
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Fig. 13. Baryon number densities for hybrid star matter in¬ 
cluding a first-order phase transition (£ = —1.29). One can 
observe the increased strangeness in the high-density phase. 


merited by the appropriate boundary conditions. These 
are given by demanding a vanishing heat flux at the cen¬ 
ter of the star and fixing the luminosity at the surface. 
Here, one has to take care of the relation between the 
temperature of the mantle and the actual photosphere as 
described in mm- The calculation takes into account all 
standard cooling channels via neutrino emission, includ¬ 
ing direct and modified Urea processes as well as neutrino 
bremsstrahlung and the pair-breaking-formation process 
that occurs during the onset of pairing. The pairing chan¬ 
nels include neutron singlet ( 1 5'o) in the crust, neutron 
triplet P 2 ) in the core as well as proton singlet ( 1 5'o) in 
the core. In addition to that one may also have to con¬ 
sider the different quark pairing that may occur in the de- 
confined quark phase. Among the most promimnnt quark 
pairing patterns are the Color-Flavor-Locked (CFL) and 
the 2-color (2SC) superconductivity. 

As an example of such a study, Fig. [14] shows the 
cooling curve of the two equal mass stars, one from the 
normal branch and one from the twin branch, which are 
shown in Fig. [12] Although both stars have equal mass of 
M = 1.68 M©, their thermal evolution is distinctly differ¬ 
ent from each other. As can be seen, the strangeness-rich 
twin star cools significantly faster than the star from the 
ordinary branch. This might help to distinguish stars that 
might have similar masses but a very different structure. 
Note, however, that this result depends on the energy gap 
used for the quark pairing. A large gap would make both 
cooling curves more similar. 



Fig. 14. Cooling curve of the hybrid star and its high¬ 
strangeness twin solution with the same mass of 1.68 Mq. As 
can be seen in the figure, the twin star cools significantly faster. 


2.5 Including baryonic parity partners 

With respect to hadronic models, there has been a number 
of studies looking at nuclear matter properties and also at 
stellar matter with a different formulation of chiral sym¬ 
metry and the baryonic fields. These refer to the so-called 
parity-doublet model [5T1I521I53] . The main feature of the 
doublet model is the inclusion of both, the positive and 
negative parity states of the baryons. By arranging them 
in doublets, the components of the fields defining the par¬ 
ity partners and c p _ transform differently with respect 
to chiral transformations: 

P+r = Rp+R P+l = l P+l , 

if'_ R = Lip-R p- L = Rp~L , (20) 

where L and R are rotations in the left- and right handed 
subspaces. This behaviour is very similar to the mesonic 
analogue of the sigma meson and the pion as grouped to¬ 
gether in the linear sigma model. In the latter case, the 
combination of sigma and pion a 2 + 7r 2 is chirally invari¬ 
ant, whereas the separate terms are not. Here, including 
the opposite parity states helps to achieve the same ef¬ 
fect. It allows for a chirally invariant mass-like term in 
the Lagrangian of the form: 

m 0 (< P - 75(^4. - ( p + - y 5 ip _) = 

'moi’P-Lf+R ~ P-rP+l - P+lP-r + P+rP-l), (21) 

where mo is a mass parameter. 

The SU(3) extension of the parity-doublet approach 
has been developed in Ref. [54]. In general one can write 
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the following SU(3)-invariant baryon-meson interaction 

C B = Tr(Si^S) + m 0 Tv((S l5 T 2 S) + D^Tr(S {S, ~ }) 
+ fPTi(S [S, »]) + S^Tr(S)Tr(SS) 

+ Df Tr(Sr 3 {V ~ }) + F^Tv(St 3 [2, S]) 

+ 5f)Tr(r)Tr(^r 3 H) + D v Tr(S lfi {V», 5}) 

+ F v Tr(^ [V», S]) + S v Tv(V»)Tr(S^~) . (22) 

Here, E denotes the baryon octet that now also contains 
the doublets of all the octet states. The r matrices act on 
the doublet spinors. £ and V^ represent the scalar and 
vector multiplets, respectively. Note that there are now 
two sets of couplings F^\D^\S^ with j=l,2 for the 
usual SU(3) invariant expressions. 

After diagonalization in the doublet-spinor space, the 
following Lagrangian results 

C B = (Bim^Bi) 

i i 

+ + g P ip M . (23) 

i 

The coupling constants for the baryons Bi with the vector 
mesons cj, p and the strange meson 0 are analogous to 
Eq. <§. The main difference to the previous approaches 
appears in the expression for the effective baryon masses 
(neglecting the isovector scalar meson, for simplicity) 



Fig. 15. Phase diagram of the SU(3) parity-doublet model in 
the plane of temperature and quark chemical potential. Two 
first-order phase transitions for isospin-symmetric matter can 
be observed, corresponding to liquid-gas and chiral transitions. 
The bands correspond to a change in value of the ^ and a fields 
between 20 and 80 percent of their maximum value. 


(aiF + g^iO 2 


(too 




(24) 


where there are two sets of coupling constants to the 
baryons i g J ai ^g J ai , which are combinations of the original 
parameters. The signs db distinguish the cases of opposite 
parity baryons. As can be seen in the expression, in the 
case of vanishing scalar fields, the various doublets are 
degenerate but not massless. The expression contains a 
SU(3) breaking mass term with m s = 150 MeV, respon¬ 
sible for the generation of an explicit mass corresponding 
to the number of strange quarks n s in each baryon. 

With these relations at hand, we can perform a calcula¬ 
tion analogous to the one in Section [274] including quarks. 
The resulting phase structure of the model for isospin- 
symmetric matter is given in Fig. [15] As an interesting 
difference to the previous discussion, the model exhibits a 
second first-order phase transition, which is driven by the 
appearance of quarks as well as the parity-doublet partner 
of the nucleons. This transition ends in a critical end-point 
around a temperature of 60 MeV. Thus, it might be very 
useful to study the implications of such a phase structure 
in neutron star merger calculations that might sample the 
conditions close to such a critical point. 

Besides the parametrization used in the calculation de¬ 
scribed above, which reproduces a high value for the com¬ 
pressibility at saturation, we define another parametriza¬ 
tion with a lower value for the compressibility. TOV solu¬ 
tions for both equations of state are shown in Fig.[l6] All 



Fig. 16. TOV solutions for hybrid stars in the parity-doublet 
approach. Two different equations of state are compared. Re¬ 
sults for vanishing temperature and T — 30 MeV are shown. 
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Fig. 17. Particle densities as function of distance from the 
center of the star for the maximum mass star of the stiff equa¬ 
tion of state of Fig.[l6]at T=0. The complex particle cocktail 
contains the octet of baryons and their parity doublets as well 
as quarks and leptons. 



age [y] 


Fig. 18. Stellar temperature as function of stellar life time 
resulting from a cooling simulation. Observational data are in¬ 
cluded as symbols. The curves show results for different star 
masses. The simulation includes quark pairing with a gap of 
10 MeV. 


curves for the star solutions show the familiar behaviour 
of a kink signalling the crossing of the first-order tran¬ 
sition point, with a maximum between 1.7 and 2 M 0 , 
respectively, for the soft and stiffer equation of state at 
zero temperature. In addition, the figure contains results 
for a fixed temperature of 30 MeV with slightly increased 
maximum masses. Note that the fixed temperature calcu¬ 
lation is not an ideal approximation for a snapshot of a 
proto-neutron star thermal evolution, where a calculation 
at fixed entropy per baryon would have been more real¬ 
istic and should be performed in future work. The tem¬ 
perature used mimics the maximum temperature effect to 
be expected in the core of the star. It corresponds to an 
entropy per baryon of about 2 in the stellar center, which 
is a reasonable early-stage value compared to numerical 
simulations [5511561157] . 

The particle composition of these stars can be quite 
complex as is presented in Fig. [IT] for the stiff equation 
of state at zero temperature. These stars contain baryons, 
quarks, as well as the opposite parity states of the nucleons 
and some hyperons. 

Finally, we look at the temperature evolution of stars 
within this model. Fig. [18] shows results for a range of 
star masses for the stiff equation of state. The calculation 
includes a gap of 10 MeV for the quark pairing in the 
core. The figure also includes a range of observed stellar 
temperatures and star ages derived from the spin-down 
of the pulsar as well as from kinematics, tracing the star 
back to the original supernova. Overall, there is a quite 
good agreement of the different measurements with the 


region of possible temperatures. However, a simultaneous 
measurement of star masses would be very helpful in con¬ 
straining neutron star models via studies of stellar thermal 
evolution. 


3 Possible conflict with lattice results 


As has been documented in the previous sections, one 
nearly universal ingredient for all hybrid star calculations, 
within the framework discussed in the article as well as in 
many other calculations, is the inclusion of a strong repul¬ 
sive force between the quarks [58ll29ll3Q[l59ll60l[611l62j . This 
force stiffens the equation of state and allows for the sur¬ 
vival of a larger quark core without immediate collapse 
to a black hole. The calculation in Section 12.31 did not 
introduce a non-zero vector coupling of the quarks. How¬ 
ever, as it was shown, this then only allows for a region 
of a mixed phase of quarks and baryons if one performs a 
Gibbs construction of the phase transition. 

When one applies the different model approaches to 
matter at vanishing baryon chemical potential, one has 
the advantage to be able to compare results with lat¬ 
tice QCD calculations. This has been used in adjusting 
the model pa rame ters of the approaches, discussed in Sec¬ 
tions [2T3JJ274] and 2J5 , and a direct comparison was shown 
in Fig. W Overall, without too much fine-tuning a rea¬ 
sonable agreement with lattice data can be achieved for 
these star models. Unfortunately, such comparison can¬ 
not easily go further, as lattice calculations at non-zero 
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Fig. 19. Taylor coefficient of the pressure C 2 as function of 
temperature at zero chemical potential for the quark hadron 
model described in Section [2~4l The solid line was obtained for 
zero interaction, the dashed line represents results for g q0J — 
gNuj/%- Again, a small interaction strength causes results that 
are inconsistent with lattice data. 


Fig. 20. Same as Fig. 


19 


but for a standard PNJL model. The 


curves correspond to different values of the repulsive quark- 
quark interaction, G v — 0, G^/2, GV 2GV where G s is the 
strength of the scalar interaction. As can be seen, even a rela¬ 
tively small interaction strength leads to strong disagreement 
with lattice data. 


chemical potential are notoriously difficult. In this case, a 
phase in the weighting factor of the path integral leads to 
very problematic numerical cancellations in the sampling 
of the integral. Some effort has been put into a reweighting 
method to calculate observables at finite chemical poten¬ 
tial. A main result was the computation of a critical end¬ 
point of a first-order transition line at T = 162 MeV and 
li b = 360 MeV with larger statistical errors. This result 
has been used in Section [2ffi| to fix the model parameters. 
However, it is still debated whether this critical point ex¬ 
ists or not. Another way of extrapolating into the regime 
of non-zero chemical potential is to do a Taylor expansion 
of the grand canonical potential and calculate coefficients 
of an expansion in /ib- The expansion of the pressure can 
be written as 


P(T , /i b ) 

2^4 


£c„ ( r)(f)’ 


n=0 


c„(T) = 


1 d n (p(T, h b )/T a ) 


n\ d(ns/T) r 


(25) 

(26) 


0 


The advantage of this approach is that the coefficients 
c n can be calculated on the lattice as expectation val¬ 
ues of operators at vanishing chemical potential. The first 
non-zero correction term of such a study is C 2 . Calculat¬ 
ing the coefficient in the quark-hadron model discussed 
in Section |2.4| and comparing it to lattice gauge simula¬ 
tions lead to the results shown in 19 [63] . The comparison 
in the figure uses somewhat older lattice data [53] . For a 
more careful study with recent lattice data, see Ref. [55] . 
Note, however, that all qualitative conclusions remain un¬ 
changed. Various theoretical curves are shown that corre¬ 
spond to different strengths of the vector interaction. As 


can be clearly seen, a non-zero vector coupling leads to 
very significant differences with respect to the lattice data. 
In order to demonstrate that this is not a model-specific 
result, an analogous calculation was done for the PNJL 
model. These results, depicted in Fig. |20| demonstrate the 
same behaviour. A reasonably large value of the vector in¬ 
teraction coupling Gy results in clear disagreement with 
the lattice simulations. Also, an analysis of the pseudo- 
critical phase transition line as function of chemical po¬ 
tential that, according to Ref. [55] . signals a strong vector 
interaction term, does not change this picture. As it has 
been shown in Ref. [55] . the rather flat transition line is 
governed by the properties of matter below the transition, 
i.e. the hadronic matter that naturally has a strong vector 
interaction, but not by the quark phase (see also Refs. m 
155]). Therefore, this problem still persists, at least for the 
models currently under investigation. 


4 Conclusions 


In this article we have discussed several approaches con¬ 
necting hadronic descriptions and quark models in a uni¬ 
fying way, which allows for first-order and higher-order 
phase transitions as well as crossovers. Such models can 
reproduce a mixed state of quarks and hadrons, even with¬ 
out the appearance of any sharp phase transition in dense 
matter. As one interesting result, we observe twin star so¬ 
lutions for some model parameters, yielding two sets of 
hybrid stars with different strangeness content and radius 
but with the same mass. As was shown, these stars would 
have a distinctly different thermal evolution. 
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In accordance with other theoretical studies, we find 
that in order to obtain a star with a large quark core or 
a more extended mixed phase of quarks and hadrons, the 
quark interactions have to be strongly repulsive. However, 
such assumption leads to problems in the agreement with 
lattice data at small densities (that study the Taylor co¬ 
efficients for an expansion of the pressure with respect to 
baryon chemical potential). This problem should be ad¬ 
dressed before one can make more definite claims about 
the feasibility of hybrid stars with large quark cores. Hope¬ 
fully, this point can be reasonably addressed without the 
need to render the strong interaction models too complex 
with the inclusion of many additional phenomenological 
terms. 
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